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I. INTRODUCTION 



The recent observation of CP violation in the 5-meson system, by dint of the measurement of 
a nonzero, CP-violating asymmetry in the decay B°(B°) — > J/ipKs and related modes [1], opens 
a new era of discovery. At issue is whether or not CP violation is realized exclusively through a 
single phase, 8km, in the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [2], as in 
the standard model (SM). The observed asymmetry, realized through the interference of B°-B° 
mixing and direct decay, is in accord with SM expectations [3], but additional, critical tests of the 
SM picture have as yet to be realized. For example, the B°(B°) — > J/ipKs result, coupled with 
the hierarchical nature of the CKM matrix [4], suggests that significant direct CP-violating effects 
ought to exist in the 5-meson system. Moreover, the determined angles of the unitarity triangle, 
whatever they may be, must be universal, and they must sum to a multiple of n [5]. In this paper, we 
study a new method of elucidating the presence of direct CP violation. The empirical observation 
of direct CP violation in the B system would falsify models in which CP violation is associated 
overwhelmingly with \AB\ = 2, rather than \AB\ = 1, processes: such models are "essentially" 
superweak [6-9] . The observation of direct CP violation in the B- meson system is needed to clarify 
the nature and origin of CP violation in nature. 

Direct CP violation (DCPV) in the S-meson system is realized if \Aj/Af \ ^ 1, where Af is the 
decay amplitude forB^f and Aj is the decay amplitude for the CP-conjugate process B — > /; it 
can be established in a variety of ways. Let us consider the possibilities. Perhaps most familiar is 
the partial-rate asymmetry, Aqp, for which A CP oc T(B — > f)—T(B — > /), with T(B — > /) oc \Af\ 2 . 
No conclusive evidence for DCPV has been found thus far, although early results in certain modes 
have been quite suggestive [10-12]. Another pathway to DCPV is realized through the comparison 
of acp(f), the CP- violating asymmetry associated with the interference of B° — B° mixing and 
direct decay, for two different self-conjugate final states /. Were |acp(/)| — acp(tpKs) nonzero for 
a final state /, one would establish the existence of DCPV, as in the manner of the e' parameter in 
the neutral-kaon system. Indeed, such a difference has been termed e' B [13]. Experimentally, there 
is no conclusive evidence, as yet, for a nonzero e' B [14, 15]. It is worth noting that e' B can vanish for 
a specific final state, such as 7r + 7r~, even if DCPV is present in nature [16, 17]. Direct CP violation 
can also be established through the study of the angular distribution of B decays into two vector 
mesons [18-22], though no empirical limits on the DCPV terms exist as yet. Finally, DCPV can be 
established through a population asymmetry in the Dalitz plot associated with the untagged decay 
rate into self-conjugate final states [23, 24], as we will develop in detail. The last two methods 
have common elements; in particular, neither method requires tagging nor time- dependence to find 
DCPV in the neutral 5-meson system. 

In this paper, we study the method proposed in Ref. [23] to search for DCPV in the decays of 
neutral, heavy mesons. The method uses untagged, neutral meson decays into self-conjugate final- 
states [25-27] containing more than two hadrons 1 . The multi-particle (n > 2) final states realized 



A note on language: a two-particle, "self-conjugate" state is a CP-eigenstate, whereas a multi-particle (n > 2) 
state is not a CP-eigenstate, in general. Nevertheless, we shall persist in using the phrase "self-conjugate" to 
describe multi-particle states which are "CP self-conjugate states in particle content" [28]. 
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in heavy-meson decays possess a rich resonance structure, entailing the possibility of detecting 
DCPV without tagging the flavor of the decaying neutral meson. For example, DCPV can occur 
if we can separate the self-conjugate final state, via the resonances which appear, into distinct, 
CP-conjugate states. This condition finds it analogue in stereochemistry: we refer to molecules 
which are non-superimposable, mirror images of each other as enantiomers [29], so that we term 
non-superimposable, CP-conjugate states as CP-enantiomers [23]. In B — > 7r + 7r~7r° decay, e.g., 
the intermediate states p + ir~ and p~7i + form CP-enantiomers, as they are distinct, CP-conjugate 
states. As a result, the untagged decay rate contains a CP-odd amplitude combination. The 
empirical presence of this CP-odd interference term in the untagged decay rate would be realized 
in the Dalitz plot as a population asymmetry, reflective of direct CP violation [23]. The use of 
untagged decays to search for direct CP violation is also possible in B — > ViV 2 decays; its practical 
advantage is that there is no loss of statistics due to tagging. 

We focus in this paper on the application of this method to untagged 5-meson decays into self- 
conjugate final states of three pseudoscalar mesons. In Sec. II, we detail, after Ref. [23], how direct 
CP violation can arise in untagged, neutral- -B-meson decays. In Sec. Ill we discuss how the expected 
size of the population asymmetry compares, on general grounds, to that of the analogous partial- 
rate asymmetry. In Sec. IV, we make explicit numerical estimates of the population asymmetries in 
particular channels and discuss strategies for their experimental realization. In specific, we consider 
the decays B d — > it + it~it , D + D~it , as well as the analogous decays B s — > K + K~tc°, DfD~ir° in 
the B s system. These final states contain the CP-enantiomers p ± ir T , D* ± D T , K*^K* , and D^DJ, 
respectively. We conclude with a summary and accompanying outlook in Sec. V. 

II. DIRECT CP VIOLATION IN UNTAGGED B DECAYS 
A. B — > (m — > 7T+7r _ 7r° 

We follow Ref. [23] and use the decay B — > pir — > it + it~it as a paradigm to illustrate how direct 
CP violation can occur in untagged, neutral heavy-meson decays. We express the amplitudes for 
B°(B°) -> p + 7r-,p-7T + ,p°n decay as 

= a kl e* p ■ p n and M pkw i = a kl e* p ■ p„ , (1) 

respectively. Under an assumption of p dominance, the amplitudes for B° (B°) — > 
n+ i.P+) 7r~(p_) tt°(Po) decay can be written as [30, 31] 

M 3w = a + _ f + + a_ + f_ + a 00 /„ = a f + a u f u + a n f n , 

(2) 

-^Stt = a+_ /+ + a_ + f_ + a 00 f = a g f g + a u f u + a n f n , 

where we have summed over the p polarization. We define 

a g = a + _ + a_ + , a u = a + _ - a_ + , a n = 2 a 00 , 

a g = a + _ + , a u = a + _ - a_ + , a n = 2 a 00 , 
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and 



/± — =t§(s+_ s to) r pvr7r (s±o) ' /o — l ! 



s_ n - s 



O °+07 fm-K 



>+-) , 



(4) 



2 /g — /+ + /- ; 2 fu — f+ f- 1 2 / n — /q , 

where we have neglected the charged-neutral pion mass difference. Note that r p7r7r (s) is the p — > 7r7r 
vertex function, which can be determined from e + e~ — > 7r + 7r~ data, and that s fci = (p k + Pif- Let 
us consider the untagged decay width: 



3vr 



ds +0 ds_ 



256tt 3 



in 



BO 



(\M 3 J + \M 3 J) =rS + r£» + r« 



(5) 



We can separate this width into pieces which contain CP-even, as well as CP-odd, amplitude 
combinations, so that 



r£ = /d$<j ^ (| a J 2 + |aJ 2 )|/J 2 + 2Re (a> n + a* g a n ) f* g f n 



(6) 



. K=g,u,n 



whereas 



F (2) 
1 3tt 



| d$ 2 Re [{a* g a u + a* g a u ) /;/„] , rg = | d$ 2 Re [«a n + «X) /„*/„] , (7) 



with 



d$ 



dg+o ds_ 
2567T 3 m^o 



(8) 



It is worth emphasizing that the are known functions, as T p7Tn (s) itself can be determined from 
data. Consequently, the coefficients of f^f l are observables: they can, in principle, be determined 
from the Dalitz plot of the 7r + 7r~7r° final state [30, 31]. 

We now consider the transformation properties of the amplitudes under CP. Under the con- 
ventions we detail in Sec. II B, we have 



CP _ 
a, < — > a 



■+ ' 



CP _ 
a , < — > a 



+- ' 



CP _ 
a oo < — > a 



00 > 



(9) 



which translate into 



CP _ 
a g < — > +a g , 



CP 

a„, < — > —a,. 



CP _ 
a n < — ► +a„ . 



(10) 



It follows that the combination a*a n + a*a n in Eq. (6) is CP even, whereas the combinations 



a* g a u + a* g a u and a* u a n + a* u a n 



;n) 



in Eq. (7) are CP odd, as first noted by Quinn and Silva [31]. The discernable presence of either of 
the last two combinations signals direct CP violation. These quantities can be measured with the 
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aid of the Dalitz plot of s +0 vs. s_ , as illustrated in Fig. 1. To see how this is useful, note that 
under p + <-> p_ we have, from Eq. (4), that 

fg,n{ S +Qi S -o) = — fg,n{ S -Qi S +o) ? /m( S +0> S ~o) = + /u( S -0> S +o) J (12) 

which implies that f*f u and f*f n are both odd under the interchange s +0 <-> s_ , whereas f*f n 
is even. Thus the CP-odd terms make no contribution whatsoever to the total, untagged decay 
rate; the behavior of Eq. (12) ensures that these contributions vanish identically in the integral over 
phase space. However, they do generate a population asymmetry about the s +0 = s_ line in the 
37r Dalitz plot. Explicitly, noting Eqs. (6) and (7), we have 

r { £[s +0 >s_ ] = +r£[ s+ o<s- ] , r£ 3) [ s+0>s _ ] = -vff[s +0 < s _ } , (13) 

where the inequalities within the square brackets indicate the regions of integration within the 
Dalitz plot. The corresponding asymmetry can be expressed as 

A _ r 37rKo> g -o] ~ r 3^[ g +0< g -o] = r^Ko^-p] + ?S [ g +0> g -o] 

3 " " r 37r [ s+ o>s_o] + r 37r [ s+0 < s _ ] " r«[ s+0 > s _ ] ^ 

A nonvanishing value of A 3n reflects direct CP violation in the B — > pit amplitudes. The individual 
asymmetries, 



,(2) = lfJ[s +0 >S_ ] (3) r ( V[s +0 >S_ V : . 

° n ~ T (1) \s >s 1 ' ° n ~ T (1) \s >s 1 ' 1 J 



are distinguishable through their location in the Dalitz plot; moreover, they are proportional to 
different combinations of weak amplitudes, as clear from Eq. (7). The fa functions embedded in 
the definitions of r^ 3 "* ensure that the population asymmetry is restricted to the p k bands. In 
specific, r 37r is nonzero if the population of the p + band for s_o > s+o is different from that of the 
p~ band for s_ < s +0 , whereas is nonzero if the population of the p + -p° interference region, 
which occurs for s_ > s+o, is different from that of the p~-p° interference region, which occurs for 
s_o < s +o- It turns out that T 3 J is sufficiently small, on general grounds, that its impact on T 3 J is 
negligible. 



We have illustrated how DCPV can occur in the untagged decay rate. The untagged decay rate 
is distinct from the untagged combination of time-dependent rates. To appreciate the consequences 
of this, we consider the time dependence of neutral, B- meson decay. We let \B°(t)) denote the 
state which is tagged as a B° meson at time t — 0. We define the mass eigenstates under the 
weak interaction as \B H ) = p\B°) — q\B°) and \B L ) = p\B°) + q\B°) ; we denote their masses 
and widths by M H L and T HL , respectively. Denoting Aif(t) (Mf(t)) as the decay amplitude for 
B°(t) (B°(t)) — > / decay, and assuming mixing CP- violation to be negligible, so that \q/p\ = 1 , 
we have [32, 33] 

Art 



M f (t)\ 2 = e- rt \ (JM/ + \M s f) cosh^- + \ (\M S \ 2 - \M s f) cos(AMt) 



+ Re (^MfM^J sinh 



ATt Im ( ^M}M f ) sin( AM t) 



(16) 
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FIG. 1: The allowed Dalitz region for B°,B° — * 7r + 7r~7r decay. The decays mediated by the p^ir^ 
intermediate states dominate the Dalitz plot; their location is indicated schematically by the shaded bands. 
The p° band, reflective of the /?°7r° intermediate state, is indicated schematically by the dotted lines. 
An asymmetry in the population of the Dalitz plot about the s+o = s -o "mirror line" signals direct 
CP violation. Note that the asymmetries A^J and A^J , as given in Eq. (15), describe the population 
asymmetry of the /9 ± bands and that of the overlapping p^-p° bands about the mirror line, respectively. 



and 



\M f (t)\ 



-rt 



Mf\ + \M f n cosh 



+ Re ( -M)M f ) sinh 



V 



Art 

ATt 



MA 2 ) cos(AMt) 



1m ( ^M}M f ) sin ( A.I/ 



(17) 



We have used the conventions of Ref. [33], so that 2 T = T H + T L and AM = M H — M L , though 
we define Ar = T H — T L . The untagged combination of time-dependent rates is given by 



\Mf(t)r + \M f (t)\ 



-rt 



\M f \ z + \M f \ z ) cosh— + 2 Re (~M* f M f ) sinh — 



V 



The population asymmetry associated with the untagged combination of time-dependent rates need 
not be a signal of direct CP violation at nonzero t. The population asymmetries are generated by 
terms of form f*f u or f*f n , as in Eqs. (7,15); at t ^ 0, A4*jAij- enters as well. Such terms are 
associated with CP-odd amplitudes in Re(Ai*j-M. *) , but with CP-even amplitudes in lm.(A4%A4f). 
The pieces of interest in A4 *tMf are of form 



M}M 3 



gu ■ 



where a^ u = a*a u ± a* u a g and the replacement g,u <-> u,n yields the /„/ n terms. We note 

± cp ±* 
az, < — > =fo 



gu > 



(19) 



(20) 
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so that the CP properties of (Re, Im)(.Mj.Mj) follow. Since q/p is complex, \m(M*fM. j) enters 
R,e(qAi*fAij-/p) as well, so that the population asymmetry associated with this t-dependent term 
can be generated by direct CP violation, but does not occur exclusively because of it. In a similar 
vein, we conclude that a population asymmetry in |.M/(£)| 2 — |*M/(t)| 2 need not be reflective of 
direct CP violation; at t — it is decidedly not. Since Ar <C T in the SM [34], though it is likely 
more significant in the B s system [35, 36] 2 , the practical impact of this observation ought to be 
minor; it nevertheless remains that the population asymmetry signals direct CP violation in the 
untagged decay rate, rather than in the untagged combination of tagged, time-dependent rates. 
However, if the latter quantity is integrated in t over a symmetric region about t = 0, the undesired 
term vanishes, and the resulting time-integrated, population asymmetry is also reflective of direct 
CP violation. 

B. B — ► TT + TT-TT 

Herewith we demonstrate that the population asymmetry in the untagged combination of time- 
independent rates in B — > 7r + 7r~7r° decay is a signature of direct CP violation; we need not confine 
ourselves to the regions of the Dalitz plot dominated by the p resonance. Indeed, our demonstration 
follows for all three-body decays to final states which are self-conjugate in particle content and whose 
particles are spinless; the final-state particles need not be related by isospin symmetry. 

Let us consider the amplitude for B — > 7r + 7r~7r° decay. We enumerate the possible contributions 
by the orbital angular momentum of the final-state mesons, as this suffices to characterize the CP 
of the final, three-particle state. We note, in passing, that a complete parametrization of the decay 
of a pseudoscalar meson to 7r + 7r~7r° in terms of the isospin I of the pions and |A7| of the weak 
transition can be found in Ref. [37]. The final state must have zero total angular momentum, so 
that we write (^1^2) 1(^3) v where / denotes the total orbital angular momentum quantum number 
of the particles in parentheses, combined to yield a state of total angular momentum zero. 3 We 
have implicitly summed over the magnetic quantum numbers, m, possible for fixed I. The decay 
mechanism can distinguish the way in which the three pions are coupled to yield a state of total 
angular momentum zero. For example, in decays mediated by a two-body intermediate state, such 
as B — > p + 7r~ — > (n + n ) 1 ( y n') 1 or B — > p~n + — > (7r _ 7r°) 1 (7r + ) 1 , the differing ways in which two 
of the three pions are coupled to a 1=1 state reflect distinct weak amplitudes. We denote the 
amplitude for B — > ( y n + 7i ) 1 (n') 1 decay as -M + _(p + ,p_,p ), so that the first subscript denotes the 
net charge of two pions coupled to a / = 1 state and the second subscript denotes the charge of the 
third pion; its arguments p + ,p-,po denote the momenta of the 7r + , 7r~, and 7r° mesons, respectively, 



2 We note AT s d M /T d w 3 x 1CT 3 , whereas ArS M /r s = 0.12 ± 0.06, as per Ref. [33] and references therein. 

3 Note that (n 3 ) l denotes the orbital angular momentum of n 3 with respect to the center-of-mass of the 71-^2 system. 
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so that we define 



eff 



<2 



5° 



A4(p+,P-,Po) = (in+ip^n-ip^y^ H cS B°) = ^ A* ((tt + 'K~) 1 7r° 



(21) 



5 U , 



with analogous quantities for £? decay 

A^U(p + ,P_,Po) = (k T (PT) 7r °(Po)] i vr ± (p ± ) ( 



n 



off 



5 ( 



0, 



5 ( 



^(p + ,p_,Pb) = ([7r + (p + )7r-(p_)],7r°(po). Kfr 5°) = E **((* + 0, Qi B °) > 



(22) 



where we write the \AB\ = 1 effective, weak Hamiltonian as 

K s = E( A A + A*Qj) • 

i 

To determine the transformation properties of these amplitudes under CP, we have 



(23) 



CP\B) = Vb \B) , C7V) 



7T 



CP\n + ) = , 



cpq^cp^ = ^Ql , 



(24) 



where Vb,tt,q are arbitrary phase-factors. Armed with these conventions, we turn to the determina- 
tion of the CP-properties of the three-meson final state. Working in the rest frame of the two pions 
coupled to angular momentum I, we note that 



P 



7T, 



Iip)^2(-P)]inip')l) = ~ Mp)^2(-P)\i^3(p')i) , 



(25) 



where the minus sign emerges from the negative intrinsic parity of the pseudoscalar mesons. Under 
charge conjugation, 



C 



7r+(p)7r-(-p)],7rVU = [n~(p) TT+f-p)], n (p% ) = ^(p)], ttV). 



(26) 



[7r+(p)7r-(-p)],7rV)*> = ~C [n^n-i-p)}^'),} = (-1)< +1 7r + (p)],7r°(p), , 

CP 



It follows that 
CP 

[7r-(p)7r°(-p)U + (A) = - \^{P)A-P)]^-{P)1) • (27) 
We thus determine 



MP-K(Po)W(p + )z ^ B°) = ([K-(P-)APo)U + (P + )i (CPyCPQ^CPyCP 



B 



= -V*bV Q ([k + (P-)Apo)]i*-(p + )i Ql B°) , 



(28) 



where we work in the frame in which p = p, p = —p, and p + = p' . Moreover, 



7T+7T-), 7T? 



(CP^CPQi (CP) 1 CP 



(29) 



where p_ = p, p + = —p, and p = p' . The determined CP properties are independent of frame, 
so that we conclude, more generally, that 4 



M l _ + (p + ,p_,p ) = -^Q^A l ([7r + (p_)7r°(p )] i 7r (p + ),Ql 

i 

M l + _( P+ ,p_, Po ) = -^Q^A,([7r-(p + )7r°(p )]^ + (p_)^|gi 

i 

M l 00 (P + ,P-,P ) = V*bV Q (-1) 1+1 ^A i ((7r + (p + )7r-(p_)) z 7r°(po) 



B ( 



B ( 



(30) 



Ql 



B l 



or that 



M l _ + (p + ,p_,p ) 



CP 



-V*bVqM 1 + _(p_,p + ,p ) , 



-M+-(p+,P-,Po) < — ► ~V*bVqM_ + (p^,p + ,p ) , 



(31) 



•A^oo(p + ,P-,P ) 



CP 



-V*bVqM 1 00 (p^,p + ,p ) , 



where the (—1)' in the last line of Eq. (30) has been absorbed through the exchange p + <-> p_. 
Thus we observe that the CP transformation is tied to the mirror transformation of the Dalitz plot, 
P+ P- - 5 Generalizing the definitions of Eq. (3), we introduce 



M l g (p + ,p_,p ) 



I \M l + _(p + ,p_,p Q ) + M l _ + (p + ,p_,p ) 

+ {-I) 1 M l + _(p_,p + , Po ) + (-l) l M l _ + (p_,p + , Po ) 
\ M l + _(p + ,p_,p )+M l _ + (p + ,p_,p ) 

- (-l) l M l + _(P-,P + ,Po) ~ (-l) l M l _ + (p_,p + , Po )] , 
M l n ( P+ ,p_, Po ) = M l 00 (p + ,p_, Po ) , 



M l u (p + ,p_,p ) 



(32) 



4 We can check the CP assignments by determining the CP of associated, two-body modes. Noting 
CP\p°) = +\p°) and CP\p+) = n p \p~), we find ( 7 : (p )p (p p )\Q l \B") = r lB r ]Q (7T (po)p ( Pp )\Ql\B ), as 
well as (n + (p + )Tr-(p-)\H st r\p (p P )) = (tt+ {p + )-K~ (p-)\nl tr \p {p p )) , where p p = p+ + P- and p p = 0. 
Moreover, (tt+ (p + )p- (p^Q^B") = -V* B VQV* P vA^~ (p+)p + (p p )\QI\B°) and (tt" (f>_)7r°(p )|Ktr|p- (>„)) = 
?7p?7;(7r + (]j_)7r (po)|Wl tr | ( o + (pp)), where p p = p_ + p and p p = 0. Of course H str = U\ tr . Removing the 
matrix elements associated with p — ► 7T7r, we see that £? — > pit transformation properties extracted from Eq. (30) 
are in accord with our direct calculation. Choosing, in addition, rj B rjQ = +1 and 77*77,,- = — 1 yields Eq. (9). 

5 We can use the amplitudes of Sec. II A to check Eq. (31). Choosing rf B r\Q = +1 and noting / ± P t — § — / 
and / p t^~ -/ Q , so that M 37I (p+,P-,Po) = a +- f+ + a_ + f_+ a m f does yield M 3v (p + ,p_,p ) =a + _f + + 
a |_/_ + a 00 / upon CP transformation, as per Eq. (2). 



9 



and 



M l g (p + ,p_,p ) 



M l u (p + ,p_,p ) 



M l + _(p + ,p^,p ) + M l _ + (p + ,p_,p ) 

+ (-l) l M l + _(p_, P+ , Po ) + (-l) l M l _ + (p_,p + , Po ) 
M l + _(p + ,p_,p ) + Mi + (p + ,p_,p ) 

- (-l) l M l + _(p_,p + ,p ) - (-l) l M l _ + (p_,p + ,p ) 



M l n {p + ,p_,p ) = M l 00 {p + ,p_,p ) , 
where we note 6 



(33) 



CP 



M l g (p + ,p_,p ) -v*bVqM 1 9 (p_,p + ,p ) , M l u (p + ,p_,p ) ^ ~V*bVqM 1 u (p_,p + ,p ) , 



CaI i \ CP 

M n {p + ,p_,p ) < — ► 



(34) 



-V*BVQM l n (p_,p + ,p ) . 
To exploit these transformation properties, we parametrize the B — > 7r + 7r~7r° amplitude as 

oo oo 
M 37T ( P+ ,p_,p ) = J2( M +-+ M - + + Ml 00) = T,( M 9 + Ml u + M n) , 



(35) 



1=0 



1=0 



where we have suppressed the common arguments throughout, so that M- l ± T is, indeed, 
■M l ± Zf (p + ,p_,p ). Analogously, for the B — > 7r + 7r~7r° amplitude we have 



(36) 



1=0 



1=0 



A two-body decay mechanism distinguishes the terms of fixed /; however, our parametrization is 
not limited by this assumption. Note that M l + _(p + ,p_,p ), e.g., includes the product a + _f + of 
Sec. II A. With these definitions, the untagged decay rate can be written as 



\M 3 J 2 + \M„\ 



E E 

K€g,u,n j,l=0 



(M{*M l K + Mi*M l K ) + 2 MMi*M l x + M j K *M l x ) 



, (37) 



where A G g,u,n and A > k means that the g,u,n subscripts are not repeated. Applying Eq. (34), 
we observe that 

\M 3n (p + ,p_, Po ) | 2 + \M 37T (p + ,p_, Po )\ 2 ^ \M 37T (p„,p + ,p )\ 2 + \M 37T (p_,p + ,p )\ 2 , (38) 

so that the untagged decay width, T 3n , realized from integrating Eq. (37) over the invariant phase 
space, as in Eq. (5), is a manifestly CP-even quantity. This follows because the integration region 



6 We can use the amplitudes of Sec. II A to check Eq. (34) as well. Choosing rj^r/Q = +1, we see f g V t — % —f g and 
fu P t^5~ +/„, so that M 3n (p +7 p_,p Q ) = a g f g + a u f u + a n f n docs yield M 3 „(p + ,p_ 7 p ) = a g f g + a u f u + a n f n 
under CP transformation, as per Eq. (2). 
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is itself p + <-> p_ symmetric. The population asymmetry, however, allows us to observe direct CP 
violation. Forming 

a _ r 37r [s +0 >s_ ] ~ r 37 i-[g + o<s_ ] 

377 ~ ^[s +0 >s^]+T 3 Js +0 <s_ ] ' 1 J 

where 

r 3>+o> s -o] - r 3 >+o< s -o] = / d$ { [\M 37T ( P+ ,p_,p )\ 2 + |-M 3u (p + ,p_,p )| 2 ] 

•/s +0 >s_ 

- \M 37T (p_, P+ ,p )\ 2 + \M 37V (p_, P+ ,p )\ 2 ]} > (40) 

and noting Eq. (38), we observe that A 37T is odd under p + <-> p_ exchange and that it is also 
manifestly CP odd. Consequently, the observation of a population asymmetry, that is, the failure 
of mirror symmetry across the Dalitz plot of the untagged decay rate, signals the presence of direct 
CP violation. To determine the amplitude combinations which give rise to a population asymmetry, 
let us consider the transformation properties of the amplitudes of Eq. (37) under p + <-> p_, as this 
is the mirror transformation of the Dalitz plot. We have 

M\ P trZ- (-i)'A^ , M l u » -(-l)'A*!, , M l n » (-l)X , (41) 

where analogous relationships hold for the barred amplitudes as well. We thus determine that if 
j + 1 is even only the amplitude combinations 

M j g *M l u + M j g *M l u and Mi*M l n + Mi*M l n , (42) 

are odd under the mirror symmetry of the Dalitz plot, whereas if j + I is odd only 

M j g *M l g + M j g *M l g , 
Ml* M l u + Ml* M l u , 

Mi* M l n + Mi* Ml , and Ml* Ml + M j g *M l n , (43) 

are odd under the mirror symmetry of the Dalitz plot. These amplitude combinations generate a 
population asymmetry in the untagged decay rate, which we have established as a signal of direct 
CP violation. As a plurality of amplitude combinations exist, they can act in concert to dilute the 
population asymmetry we have discussed. Nevertheless, the observation of a population asymmetry 
in the untagged decay rate to three pseudoscalar mesons is an unambiguous signal of direct CP 
violation. 

Our discussion generalizes the CP-odd amplitudes found under an assumption of p dominance, 
Eq. (11), in two ways. Firstly, Eq. (42) shows that combinations analogous in structure to Eq. (11) 
exist for all even j + 1. Equation (43) shows that direct CP violation can also be generated for odd 
j + I, such as through s-wave and p-wave interference. 
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III. COMPARISON TO THE PARTIAL-RATE ASYMMETRY 



We wish to compare the size of the population asymmetry in the untagged decay rate with that of 
the partial rate asymmetry in the same decay channel. We will show that the two asymmetries can be 
comparable in size, favoring the study of direct CP- violation in the untagged decay rate on statistical 
grounds. We begin by considering the amplitude combination a*a u + a*a u . We can write dki of 
Eq. (1) in terms of its tree and penguin contributions, as per exp(i/3) a kl = T kl exp(— ia) +P ki , where 
both T kl and P kl are complex, and the weak phases (3 and 7 are given by exp(— i(3) = V tb V t * d /\ V tb V t * d \ 
and exp(i7) = V* b V ud /\V* b V ud \ . We take V kl to be an element of the CKM matrix, and we employ 
CKM unitarity to write a + (5 + 7 = n . Noting exp(— i/3) a kl = T kl exp(ia) + P kl , we have 

* 9 = T g e--™ + P g , fl u = T t( e- to + P tt , 

a = T e iQ + P a = —T e ia — P 

9 g ' g 1 u 1 u c 1 u 1 

where T g = T + _ + T_ + and T u = T + _ — T_ + , with similar relations for P gu - We omit the overall 
factors of exp(±i/5) which ought to appear in Eq. (44) as they contribute to neither |.M/| 2 , |-M/| 2 , 
nor qM* f M f /p — we note that q/p = exp(— 2i/5) in the SM if \q/p\ = I. With r K exp(i5 K ) = P K /T K 
for k = g,u, so that r K = \PjT K \, we find [23] 



a g a* u + a g a* u = -2 T g T* sin a 



r g sin 5 g + r u sin 8 U - i (r g cos 5 g - r u cos 6 J 



(45) 



The decay width associated with this amplitude combination, Tf^ in Eq. (7), contains 

Re[(a g a: + a g a* u )f 9 f:] , (46) 

so that it is sensitive to both its real and imaginary parts, because the function f f* is also complex 
- its imaginary part is generated by the nonzero p-resonance width. If we assume T g T* to be real, 
for definiteness, and consider the imaginary part of Eq. (45), which accompanies lm(f g f*) , we see 

that rf^ can be nonzero, so that direct CP violation can exist, even if the strong phases of were 
to vanish. Merely r g or r u must be nonzero, which is easily satisfied. Consequently the observables 
of Eq. (15) can be nonzero irrespective of the strong phases of a; if the width of the resonance in 
the CP-enantiomer is nonzero. Note that were T g T* complex, the p width would not be needed to 
generate DCPV. Similar considerations exist in B — > V1V2 decays. There, too, direct CP violation 
can occur in the untagged decay rate; moreover, its presence does not demand the strong phase of 
the weak transition amplitude to be nonzero [18]. 

In contrast, the partial-rate asymmetry associated with B° — > p + n~ decay is 

v b^p+k- - r S o^ p - w + oc \a + _\ 2 - \a_ + \ 2 = -4\T + _\ 2 r + _ sin5 + _ sina , (47) 

so that both a and 5 + _ must be nonzero to yield a nonzero partial rate asymmetry. Generally 
one expects l-Mp+^-l > lA^p-^+l [38-41]; this is borne out by experiment [42]. In the limit where 
M-p-tt+ becomes negligibly small, we see that Eq. (45) becomes 

1 1 2 

a g a* u + a g a* u — > — 4|T + _| r + _ sin5 + _ sin a , (48) 
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so that the population asymmetry of Eq. (15) is identical to the partial rate asymmetry in this 
limit. 

In the following section, we illustrate these ideas numerically and discuss their experimental 
realization in greater detail. After our analysis of B — > pit — > 7r + 7r~7r° decay we turn to other modes, 
studying B -> D* ± D T -> D+D-jt , B s -> K*^K^ -> K+R-Tt , and 5 S -> D^DJ -> DfD~n° 
decays. We discuss the utility of these modes in realizing tests of the SM. 



IV. NUMERICAL ESTIMATES 



In this section we estimate the population asymmetries associated with the untagged decay rate 
in a variety of decays to self-conjugate final states. We do this for definiteness, though our essential 
conclusions do not rely on these estimates: the population asymmetry in the untagged decay rate 
can be of the same numerical size as the partial rate asymmetry in the comparable, tagged decay, 
making the search for DCPV in the untagged process advantageous. 

We begin our numerical analysis by considering the effective, weak Hamiltonian for the decay 
b — > qq'q' and its CP-conjugate, where q = d,s and q' = u, d, s, c. In the SM, we have [43] 



G 



n cS = — 



F 

y/2 
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Af> (C^ + C 2 QX) + AM (C.Ql + C 2 Q C 2 ) - ^ C i Q i + C 8G Q 8G 



+ H.c. , (49) 



where G F is the Fermi coupling constant and the factor 
element Va 



V q , b V*, q contains a CKM matrix 



■ tJ . The Wilson coefficients C i are evaluated at the renormalization scale /i, and Q 1 ... 1Q 
are four-quark operators, whereas Q 8G is the chromomagnetic penguin operator. The C i and Q i 
are detailed in Ref. [43], though we interchange C l Q\ <-> C 2 Q\ , so that C x ~ 1 and C 1 > |C 2 | 
C , 3,...,io,8G • For the CKM matrix elements, we adopt the Wolfenstein parametrization [4], 



V. 



ud 



1 - A 2 /2 , 



V 



A , V ub = AX 3 (p - irj) , 



v cd = -A, = l-A 2 /2, 

V^, = AA 3 (l-p-i77) , V ts = -AX 2 



V cb = AX 2 



(50) 



V tb = 1 



neglecting C(A 4 ), where the parameters are given by [38] 

A = 0.2205 , A = 0.81 , p = 0.12 , r] = 0.34 . 



(51) 



The three-body decays we consider proceed from hadronic two-body weak transitions of the form 
B — > XY, for which the amplitude is given by 

AWy = (XY\n cS \B) . (52) 

For definiteness and simplicity, we employ the generalized factorization approximation [38, 44], 
which makes use of next-to-leading-order perturbative-QCD to estimate the needed strong 
phases. In this approach, the matrix elements C t {p) (XYlQ^pi^B) are approximated by 
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Cf Q ' {XY\Q,j\B) tree , where the Cf s are effective Wilson coefficients and the (W|<5j|-B) tree are cal- 
culated under the factorization assumption. The Cf ff enter the decay amplitudes through the 
combinations 

^eff rieE 

/-/eff | °2i n ^-/eff | ^2j-l , K o\ 

a 2i-l — L/ 2i-l + -J-T > a 2i — L/ 2i Tj— > 1 <JJ 

-''c c 

where i = 1, • • • , 5 and iV c = 3 . We adopt the a i values at \x = 2.5 GeV calculated in Ref. [38], 
collected in their Tables VI and VII. For reference, we report these numbers in Table III in the 
Appendix. Our purpose is to gain an impression of the relative size of the partial rate and population 
asymmetries in a variety of modes. 



A. B -> TT+TT-TT 

We now evaluate the CP-odd observables in untagged B°, B° — > 7r + 7r~7r° decays. To begin, we 
assume, as in Sec. II, that B°, B° 7r + 7r~7r° decays are dominated by pn intermediate states. The 
B — > pit amplitudes we employ, as well as our ancillary input parameters and form factors, are 
given in the Appendix. The population asymmetries of Eq. (15) are computed using Eqs. (6,7). We 
begin by computing 

r« = /dS /*/, , (54) 

to find 

rys +0 >s_ ] = -r 9U [s +0 <s_ ] = (-54.57- 3.05 i) x 10 10 T B0 , 

(55) 

r un K 0>s _ ] = -r un [ s+0 <s_ ] = (2.oi - 0.08 i) x io 10 r BO , 

where the arguments in brackets indicate the regions in phase space over which the integral is 
calculated. The bracketed regions can also be selected by fixing the sign of cos#, where 9 is the 
helicity angle. We define 9 in B(p B ) — > n + (p + )n~ (p-)ir° (p ) decay via 

cos 9 = P ;' P ^ , (56) 
\P-\ \Po\ 

where the primed variables refer to the momenta in the rest frame of the 7r + (p' + )7r _ (p / _) pair, so 
that p' + +p'_ = 0. Note that s_ > s +0 corresponds to cos# > 0. The Im parts of Y gu and T un are 
relatively small, so that short- distance strong phases are necessary, from a practical viewpoint, to 
observe a population asymmetry. Employing the model of Ref. [38] and the parameters reported in 
the Appendix, we find 

a* g a u + a* g a u = (0.616 - 1.683 i) x 10" 18 , a* u a n + a* u a n = (0.186 - 0.664 i) x 10~ 18 . (57) 

It is useful to report the separate contributions to the asymmetries, as per Eqs. (6,7), as well. 
To this end, we define 

r B o r B o 
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so that Eq. (15) can be rewritten 



(59) 



We have 

B£>[ S+0 >s_ ] = -7.8 x 1(T 7 , B^[s +0 >s_ ] = +6.4 x 1(T 9 , 

(60) 

B 3n = 2B£[s +0 >s_ ] = 46xl0" 6 , 

to yield 

4 2 J = -3.4% , 45 = 0-03% . (61) 

The CP- violating asymmetry of Eq. (14) is thus 

A 3w = -3.3% , (62) 

(2) (2) (3) 

which is clearly dominated by A 3lT , the T 3 J contribution. This is expected because T 3 J contains, 
in comparison, the B(B) — > p°7r° amplitude, which is "color-suppressed" at tree-level; it is also 
suppressed by phase space, as per Eq. (55). 

In constructing the population asymmetry, we need not integrate over the entire, allowed region. 
We can also study the population asymmetry over regions of the Dalitz plot; this is realized in 
Table I. Selecting the regions where the p bands overlap enhances the population asymmetry, but 
to the detriment of its detectability, as the branching ratio into the selected regions of the Dalitz 
plot, B 3n , becomes considerably smaller. Indeed, for a direct CP search it would be more efficacious 
to use the entire allowed region. 



TABLE I: The ratios = B^[s +0 >s_ ] = T%l[s +Q >s_ Q ]/F B0 for i = 2,3, total branching-ratio B 3lT = 
T 3n /T B0 , and population asymmetry A 3n , corresponding to various integration regions within the Dalitz 
plot, with m kl = sfsTi- The B 3n entries includes contributions from both sides of the s +0 = s_ line. 



Region within Dalitz plot 




fj(3) 




2B& +2B^ 
A 3n = (%) 


entire allowed region 


-7.8 x 10- 7 


+6.4 x 1(T 9 


46 x 10~ 6 


-3.3 


m P - 3r P < m +oi m ~o < m P + sr p 

m P ~ 2r P < ™+o> m -o <m p + 2r p 
m p - T p < m +Q , m_ Q < m p + T p 


-6.6 x 10~ 8 
-2.8 x 10~ 8 
-1.4 x 10~ 9 


-4 x i(r 12 

-1 x 10~ 12 

-4 x i(r 14 


22 x 10~ 7 
73 x 10~ 8 
31 x 10~ 9 


-6.1 
-7.6 
-9.3 


m p — 3F p < m H , m_ < m p + 3r p 

m p — 2F p < m H , m_ < m p + 2r p 

m p -T p < m + _,m_ <m p + T p 


-3.8 x 1(T 8 
-1.5 x 1(T 8 
-1.0 x 10~ 9 


+0.4 x 10~ 8 
+0.2 x 10~ 8 
+0.3 x 10~ 9 


27 x 10~ 7 
10 x 10~ 7 
76 x 10~ 9 


-2.6 
-2.3 
-1.9 



More realistically, B — > 3ir decays will contain contributions from other than B — > pV- 7 — > 3tt. 

(2 3) 

We can isolate the individual T 3 ^ contributions by restricting the integration regions to the p 
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bands, which lie close to the edges of the Dalitz plot, as shown in Fig. 1. To realize this, we 
integrate over the region of 37r phase-space satisfying the requirement that the invariant mass of 
two of the three pions yield the p mass within an interval of M p ± 5. This amounts to calculating 
the effective width 



rf_ p(pi+p2Mp3) = F(B -> 7r(pi)7r(p 2 )vr(p3)) 



(63) 



In an experimental analysis the strength associated with the complete p line shape would be assessed; 
our effective width Ansatz is a simple simulacrum of this more sophisticated analysis. In order to 
avoid double counting in determining the population asymmetries, we further require that the p~n + 
(p + 7r~) contributions be restricted to the s +0 > s_ (s +0 < s_ ) region. The resulting effective 
branching ratios for 5 = 0.3 GeV are 



and 



Bfo^ p - n+ [s + o>S-o\ = 4.9 x lO" 6 , £$L p+7r -[s +0 <s_ ] = 17.6 x 1(T 6 
^fo_ >p - 7r+ Ko>s_ ] = 16.2 x 1(T 6 , B^ p+n _[s +0 <s_ ] = 4.9 x 1(T 6 



tffo^oKo^-o] = 0.18 x lO" 6 , Bf ^ p0n0 [s +0 <s_ ] = 0.56 x lO" 6 
&g ^ p0 Js +Q >s_ } = 0.48 x lO" 6 , ^ p0 Js +0 <s_ } = 0.16 x lO" 6 



(64) 



(65) 



For reference, we report the two-body branching ratios and the corresponding rate asymmetries of 
the decay modes considered in this paper in Table II. We note that the combined branching ratios 
of the two-body B — > p ± n T branching ratios we report are slightly high with respect to recent 
measurements [45-47], but the estimates would be consistent with the data if the uncertainties 
in the input parameters were taken into account. Returning to our discussion, we find that the 
effective p°7r° branching ratios are much enhanced with respect to the corresponding two-body 
rates in Table II — this is caused by the large interference with the p ± 7r =F contributions. The 
population asymmetries about the s +0 = s_ line are then 

j+- _ r j^B^p-7r+[ g +0 >g -o] ~ r ff^p+ ,r-[ S +0 <S -o] = _ g ^ 
r e B B-, p - w + [ S +0 >S -o] + ^°B,B^p+7r- [ S +0 <S -o] 

(66) 

^00 _ ^B,B~> P Ojt^ S +0 >S -o\ ~ r B,B^p"7rP[ g +0 <g -o] _ _^ 
^ C B,B^p°n°[ S +0 >S -o\ + ^B^p0 7r ()[ s +0 <S -o] 

to be compared with A 3 J ,A 3 J, respectively, in Eq. (61). There is very little T 3 J contribution in 
A 3 ~, whereas in A 3 ® the dominance is much less pronounced, yielding ~— 2.6 in the —4.6. 
Practically, it is challenging to measure A®®, due to the smallness of the effective //V° rates. We 
note that the numerical value of A 3 ~ is fairly close to that of A CP for B° — > p + 7r~ in Table II, 
and this arises from the dominance of the B° — > p + 7i~ and B° — > p~7i + contributions in the p^ 
regions of the Dalitz plot. 

We expect that B, B — > 7r + 7r~7r° decay is not populated exclusively by B,B — > pit contribu- 
tions; other resonant and nonresonant contributions ought to occur — and can potentially enter 



16 



TABLE II: Estimated two-body branching ratios and rate asymmetries using the model of Ref. [38] and 
the input parameters given in the Appendix. 



Decay mode _ T b^vp ~ t b^vp 

B^VP ° Bb - vp ° B *^P ^cp-y——TY^^ 



B°- 




23.2 


21.4 


-4.1 


B a - 




6.47 


6.30 


-1.3 


B°- 


->p°7r° 


0.047 


0.068 


19 


B° — > 


D*+D- 


376 


384 


1.0 


B° — > 


D*-D+ 


322 


323 


0.3 




K*+K~ 


3.44 


4.76 


16 




K*-K+ 


0.271 


0.375 


16 




D*+D 7 


10044 


10028 


-0.08 




D* S ~D+ 


8628 


8623 


-0.03 



the pir phase space as well. Among the latter, contributions mediated by the B* meson may be 
important [49, 50]. Such effects are expected, however, to be suppressed by the BB*ir vertex, as the 
B* is highly off-mass-shell in the relevant kinematical region [51], so that we shall not consider these 
effects further. As for contributions from other resonances, we expect the "a," or / (400 — 1200), a 
broad isospin 1 = and spin J = enhancement in tttt scattering, to play a modest role [50-52]. 
The peak of this enhancement is close to the p in mass, and the crn intermediate state contributes 
preferentially to the p°ir° phase space. To explore how corrections to p "dominance" can impact 
the population asymmetry, we assume that the a resonance is the only additional contribution to 
the 3n final states, in the p bands on the Dalitz plot. 

We express the B — > an amplitudes as 



M, 



(67) 



In the presence of the a, the B — > 3n amplitudes in Eq. (2) become 

M 3n = a f + a u f u + a n f n + a a f a , M 3n = a f + a u f u + a n f n + a a f a , (68) 



where j a is the form factor describing the a — > n + n transition. Expressions for the B 
amplitudes and j a are given in the Appendix. The resulting B — > 3n width is 



an 



3tt 



= /»f E (l«JV| 5 j^)|/j^2Re[( aX + a ; s „)/;/,, + ( a x + 5 X)/„x 

^ K=g,u,n,a 



, F (2) , r (3) p(4) p(5) 

37T ' 37T ' 37T ' 37T 1 

where are given in (7) and 



^(4) 
3tt 



= J d$ 2 Re [(a* g a a + a* g a a ) f* g f a ] , r£ = J d$ 2 Re [(a* n a a + a* n a a ) f*J a ] 



(69) 



(70) 
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With the inclusion of the a, we encounter more CP-violating observables. Under CP transfor- 
mations we also have 7 

<>><, " ~K , (71) 

so that a* u a a + a* u a a in Eq. (69) is CP even, whereas 

a* g a a + a* g a a and a* n a a + a* n a a (72) 

in Eq. (70) are CP odd, the latter serving as specific examples of the CP-violating amplitude 
combinations given in Eq. (43). As we have discussed in conjunction with Eq. (11), a nonzero value 
of any one of the last two combinations signals direct CP violation. Noting that f a is even under 
the interchange s +0 <-> s_ , we can see that both f*f a and f*f a are odd under such an interchange. 
It follows that 

rt 5) Ko> 5 _ ] = -rt 5) [ 5+0 < 5 - ] , (73) 

so that rfjt' 5 '' contribute to the population asymmetry on the Dalitz plot about the s +0 = s_ 
line, reflecting direct CP violation, as we have discussed in Sec. II. We remark that the CP-odd 
observables in Eq. (72) share with those in Eq. (11) the unusual property that they do not rely 
on the nonvanishing of the strong phases in Oj. However, the combinations in Eq. (72) arise not 
from the interference of a pair of CP-conjugate states yielding the ones in Eq. (11), but from the 
interference of the p°7r° and <77r° intermediate states which have different relative angular-momenta, 
/ = 1 and 0, respectively. 

To find the asymmetries in the presence of the a, we recompute the effective branching ratios 
for 5 = 0.3 GeV, as reported in Eqs. (64) and (65). The results are 



BgL^a+o^-ol = 4-8 x lO" 6 , Bf ^ p+7T .[s +0 <s_ ] = 17.6 x 10" 6 , 



->P - 7T"' 

BgL^+Ko^-ol = 16 - 2 x 10- 6 , B^ p+7T .[s +0 <s_ ] = 4.8 x lO" 6 , 

Sfo^oKo^-o] = 0.18 x 10- 6 , tffo^oKo^-o] = 0.58 x lO" 6 , 

B%^ p0 Js +0 >s_ ] = 0.51 x lO" 6 , Bg, „ [s +0 <s_ ] = 0.16 x 10" 6 , 



(74) 



(75) 



+/£> U 7T U L -TU — UJ ' ftv^pj 

which indicate that the effects of the a in the p ± 7i T bands are really very small. The population 
asymmetries about the s +0 = s„ line are then 



l 3?r 



^ C B,B-,p-TT+ [ g +0 >g -o] ^ C B,B-+p+n- [ S +0 <g -o] 
^ e B,B-^p-7r+[ S +0 >S -o\ + r ^B_ >p + 7r -[ S +0 <S -o] 



Id d . _n_n \s i ri' > S_nl Td d , _n_n \s i n < ^S_ 



(76) 



jOO _ B,B~>p°w» l*+0^*-0\ 1 B,B^pOnO L J +0^ °-0i 
^ e B,B^p°n°[ S +0 >S -o\ + ^ C B,B-*p Q iT»^ S +0 <S -o\ 



7 Noting CP\a) = +\a), we find (7r°(p )cr(p CT )|<3*|B ) = -V* B VQ(^ (PoMPa)\Ql\B ) where p a = p+ + p- and 
p a = 0. With rjgr/Q = +1, we recover Eq. (71). 
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Comparing the A%~ values in Eqs. (66) and (76) with A 3n in Eq. (62), we see that A 3 ~ captures 
the CP- violating effects occurring in the dominant p ± 7r =F modes, as encoded in T^. In contrast, the 
value in Eq. (76) receives the respective contributions (—2.6 — 2.0 + 1.9 — 0.3)% from p(^+ 3+4 + 5 ) ; 
so that the presence of the a leads to a sizable reduction of the asymmetry. The interpretation of 
this latter asymmetry is thus rendered more difficult. 

The untagged asymmetry, in specific, A%~, has been investigated by BaBar and Belle. BaBar 
reports [46] 

peff _ peff _ 

A P cp = r !f^ P+n ~ - r f/"^ + = -0-18 ± 0.08 ± 0.03 , (77) 

B,B-^p-n+ ~"~ B,B^p+tt- 

with the recent update A P q P = —0.11 ± 0.06 ± 0.03 [48]. The s +0 < s -o criterion we impose in 
Eq. (66) has apparently not been effected in Aq P . Thus this quantity may treat events in the 
region where the /r* 1 bands overlap inappropriately. In the theoretical analysis we report in Table I, 
the impact of the overlap region, which contains an intrinsically larger asymmetry, is diluted by its 
small contribution to the total B — > p^if^ decay rate. Thus the s + o<S-o correction may not be 
essential; A P q P may well be tantamount to —A\~. The Belle measurement, in contrast, removes 
events with an ambiguous p charge assignment from their data set [47]; this is a practical realization 
of our s+o<s_o criterion. Belle reports [47] 

^=-0.38l°i?(stat)^(syst) , (78) 

where A = —Af~. 

The central values of the BaBar and Belle measurements, albeit their still-sizeable errors, point 
to much larger CP-violating effects than what we have calculated. Interestingly, more sophisticated 
calculations show similar trends. For example, A%~ has been calculated in the QCD factorization 
approach to yield the "default results" [41] 

a+- _ _ n m +0.00+0.01+0.00+0.10 / 7Q \ 
•^vr _ u - Ui -o.oo-o.oi-o.oo-o.io ■> 

where larger, more positive asymmetries, outside the stated error ranges, can occur in various 
excursions from their default assumptions [41]. Here there is no subtlety concerning the "overlap 
regions" as their assay assumes the p meson has zero-width. Using the results of Table II, we find in 
comparison that A%~ = —0.028 in this limit, so that our numerical results are of comparable size. 
Moreover, the size of the partial-rate asymmetry in the B — > p + 7r~ mode is crudely commensurate 
with the size of A 3 ~ [41], as we have found in our calculation, so that one would expect that a 
search for DCPV in these modes through the population asymmetry would be statistically more 
efficacious. 

This concludes our discussion of B — > 7r + 7r~7r° decay. In what follows we consider how the 
study of the population asymmetry in other modes can complement tests of the SM and searches 
for emergent, new phenomena. For example, the time- dependent asymmetry associated with the 
decay of a state tagged as B° or B° into a CP eigenstate / is characterized by two parameters: C f 
and S'y. The first term is generated through direct CP violation, whereas the second is generated 
through the interference of B-B mixing and direct decay. In the SM, Sj measures — rjf sin(2/3), where 
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i] = +1 (— 1) for a CP-even (odd) final state 8 , in modes driven by quark- level transitions such as 
b — > ccs , b — > sss , and b — > ccrf , modulo small corrections [53]. The determination of sin(2/3) from 
Sj/^Ks serves as a benchmark against which Sj, from other modes, can be measured. The modes 
B — > 4>K S , B — > K S K + K~ decays, mediated via b — > sss, £> — > D* + D~ decay, mediated via 
6 — > ccd, and £> s — > D* + D~ decay, mediated via 6 — > ccs, serve as examples. Of these, only <pK s 
is a CP-eigenstate. One can, nevertheless, determine Sf from the three-body decays, of which the 
vector-pseudoscalar modes are part; in these cases the coefficient of the sin(AMt) term must be 
corrected for a non- CP- violating dilution factor to yield Sf [28, 54-56]. We focus on the three-body 
decays as the population asymmetry therein allows us to test the extent to which the corrections to 
the SM tests suggested in Ref. [53] really are small. Alternatively, such observables may well help 
us establish the character of the new physics which induce deviations from the SM predictions. In 
the following we consider B — ► D* + D~ — > D + D~tt° and B s — > D* + D~ — > D+D~tt° decays as an 
explicit realization of these ideas; we also include B s — > K* + K~ — > K + K~n° decays as well to 
show that the population asymmetry can also be quite large in the SM. 

Before turning to these examples, we note that the population asymmetry may well serve as a 
useful tool in the study of B — > K S K + K~ decay. At issue is the extent to which the K S K + K~ 
final state is CP-even once the regions associated with the resonance are explicitly removed 
from consideration [57-59]. We note as per Refs. [58, 59] that the contributions of CP-even and 
CP-odd final states to B — > [K + K~] l (K s ) l decay differ in their angular-momentum character; 
their interference can give rise to a direct CP-violating population asymmetry as described in 
Eq. (43). The population asymmetry serves as an independent empirical assay of the reliability of 
the underlying assumptions of their analysis [57-59]. 



B. B -> D+D-7V 

We write the amplitudes for B°, B° -> D* + D-, D*~ D + as 

■M- B°~^D* k D l = a kl £ *D* ' Pd ) ■^■B ( >^D* k D l = ®kl £ *D* ' Pd > (80) 

where the expressions for a kl and a kl are given in the Appendix. The amplitudes for B°,B° — > 
D + {p + ) D~{p_) 7r°(p ) are then 



M DDk = a +- f+ + a -+ f- = a gfg + a u fu , 
^DDn = «+- /+ + «-+ /- = % fg + K fu ■ 



(81) 



Here, as in the 3n case, 



a g = a A + a h , a u = a + _ — a_ + , a g = a + _ + a_ + , a u = a H — a h , 



2f g = /+ + /_, 2/„ = /+-/_, 



(82) 



8 We define the asymmetry as A(t) = {T(B° -» /) - T(B° -» f))/(T(B° -» /) + -» /)) = -C/ cos(AMt) 

S/sin(AMi). 
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in this case, however, the f± are given by Eq. (110) in the Appendix. The B°,B° — > D + D it 
amplitudes lead to the decay width 



DDtt 



d$ M 



l DDn\ 



M 



DDttI 



F (l) _j_ F (2) 
DDtt ' DDtt 



where 



DDtt 



K=g,u 



It follows that 

^DDttI S +0 >S -o\ = +r^£)7r[ S +0 <S -o] > 



^D5vr[ S +0 >S -o] — ^ DDnl S +0 <S -o\ ■ 



(0 



(83) 



(84) 



(85) 



The CP-odd combination a*a u + a*a u in this case is proportional to sin (3, as can be inferred from 
the CKM factors in the B -> D* ± D T amplitudes given in the Appendix. 

To evaluate T^^, we consider the Dalitz plot of s +0 vs. s_ for the D + D~it° final states. 
Integrating over the two regions divided by the s +0 = s_ line, we find 



®dd-A s +o >s -o] 



1.0 x 10~ 6 , 



^DDn ~ 2^5 7r [s+0> S -o] 



4.8 x 10" 



(86) 



where = /T B0 for i = 1,2. The resulting CP-violating asymmetry is 



DDtt 



A 



DDtt 



^DDnl S +0 >S -o\ ^ DDn\. S +0 <S -o\ 

T ddJ s +o>*- ] + T ddJ s +o<S-o] rg^[ S+0 >s_ ] 



+0.4% 



(87) 



More realistically, the D + D~ir° final states also receive contributions from heavier charmed reso- 
nances such as the D and D 2 mesons [60]. We can isolate the D* effects, however, by limiting the 
integration regions to the bands. Therefore, we calculate the effective rates 



peff 

1 B^D*( Pl +p 2 )D(p 3 ) 



T(B^7v ( Pl )D(p 2 )D(p 3 )) 



(M D »-«5) 2 <s 12 <(M D »+«5) 2 



in analogy to the effective B — > 3ir rates in Eq. (63). Since the D* width is very narrow, we can 
choose 5 = 5MeV ~ 50r D » + , which allows us to capture the D* contributions and avoid those of 
D 02 , which are about 400 MeV heavier than D*. Thus we obtain 



Bf ^ D *- D+ [s +0 >s_ ] = 0.98 x 10~ 4 , B%^ D , +D -[s +0 <s_ ] = 1.14 x 10" 4 , 



•eff 



eeff 



(89) 



0.98 x 10" 



+ d*-d+[ s +o >s -o] — 1.16x10 4 , Bjjo^ D *+ D - [ s +o <s -o] 

where we note that our effective DDtt rates are about 3 times smaller than the corresponding 
two-body rates given in Table II, reflecting the empirical fact that ^ D *+^ DlT — 3r D , + ^ D+jr0 . 

The population asymmetry about the s +0 = s_ Q line becomes 

[s +0 >S_ ] — b^D*+D~ [ S +0< S -o] 



A 



peff _ 



DDtt 



^cff 



B,B^D*~D+ L°+0 



Sj.n>S_ 



+ P 



eff 



B,B^D*+D- i°+0 



S 4- n<S_ 



= +0.4% , 



(90) 
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which is the same as the result in Eq. (87). Since the D* is so narrow, the population asymmetry 
can also be calculated directly from the the two-body rates in Table II, namely, 

A+- _ r BQ-P— £>+ + T W^D*-D+ ~ T B^D"+D- ~ T W^D*+D- _ , n aW / q1 x 

•^ddtt — r i r i r + r " r /c ■ ^ yi J 

1 B°->-£)*-D+ ~ 1 B°^D*-D+ ~ r B°->D*+D~ ^ L B"^D*+D- 

The D*^D^ rates in this table are roughly consistent with data [61, 62]; in particular, we note the 
most precise measurement: B(B — > D* ± D =F ) = (8.8 ± 1.0 ± 1.3) • 10~ 4 [62]. In this reference the 
time-integrated asymmetry is reported as well [62]: 

A = Z°* +D ~ - Z D *~° + = -0.03 ± 0.11 ± 0.05 . (92) 
N D . +D -+N D .- D+ v ; 

The Dalitz plot associated with the decay pathway B — > D*^D T — > D + D~tt° gives rise to a CP- 
violating population asymmetry, though T D *+^ D o w + ~ 2 r o , + _ iD+jr( ,. In this case, however, as we 
have seen, the D* is so narrow that the population asymmetry is numerically indistinguishable 
from the untagged asymmetry computed from the two-body decay rates — there is no overlap 
region to treat. Consequently, the result of Eq. (92) can be compared with our estimate of the 
population asymmetry, to which it agrees, within errors. Sd*d should yield sin(2/3) to a good 
approximation [53]. 



C. B -> D+D-tv 



s 



We now consider the decays B° s , B® — > D^DJ — > D+D~ir°, in furtherance of the SM tests of 
Ref. [53]. As well known, the pertinent penguin contribution in B®, B® — > D^DJ decay possesses a 
subdominant weak phase, so that we expect its impact to be very small — Sd*d s should yield sin(2/3) 
to an excellent approximation. Were Sd*d s to differ significantly from sin(2/3), the associated direct 
CP-violating observables would yield helpful insight into the nature of the emergent phenomena. 
The requisite formulas can be obtained from those in the B — > D* ± D =F — > D ± D t -k q case by making 
the replacements 

d -> s , B -> B s , D { * } -> (93) 

and by using the b — > s and 6 — > s entries in Table III, so that we do not report them explicitly. The 
CKM factors in the B® — »■ D^DJ amplitudes are V^V^ ~ AA 2 and V^V ts ~ — AA 2 , which are 
largely real, resulting in small CP violation. The D* width has not been measured, but it is expected 
to be extremely small [63], of the order of 0.2 keV. As a consequence, the population asymmetry 
about the s +0 = s_ line on the Dalitz plot for the D+D~ir° final states can be calculated by using 
the B s -> D*D S rates listed in Table II. Thus, 

D s D s ir p iT 1 _|_F > V / 

B S ^D*-D+ ~ L B S ^D*-D+ B S ^D*+D- ~ B S —*D*+ D~ 

which is very small, in accord with SM expectations. 
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D. R -» K+K-tt 



s 



As a final example, we consider B®,Bg — > i^* ± iT T — > K + K~tt ; the pertinent formulas are 
completely analogous to those in the B — > D + D~tt° case. We need only make the replacements 
c — > -u and d — > s for the quarks, and B ^ B s and -D — > If for the mesons. The CP- 
odd combination a*a u + d*d u in this case is proportional to sin 7, as can be inferred from the 
B s — > K* ± K T amplitudes in the Appendix. 

Integrating over the two regions on the Dalitz plot for the K + K~-rr° final states, assuming no 
other contributions, we find 

^LKo> s -o] = 2.1 x 10- 7 , B KRn = 2B% 7T [s +0 >s_ ] = 2.9 x lO" 6 , (95) 
which leads to 

*k*. = r #T ± ^ ] = +15% ■ (96) 

As we have discussed, it is better to consider effective branching ratios instead. Thus, taking 
5 = ST K ,, we obtain 

Bfo_> K *- K+ [s +0 >S- ] = 0-09 x 10" 6 , B c i^ K , +K -[s +0 <s_ ] = 1.00 x 10~ 6 , 

(97) 

&m^ K --K+[ s +o> s -o] = L4 ° x 10_6 > ^ K , +K .[s +0 <s_ ] = 0.10 x 10- 6 , 

where the effective KKir rates are roughly 3 times smaller than the corresponding two-body rates 
in Table II, reflecting the isospin relation ^ K *+_ tK „ = 3 r^.+^^+^o- We thus find the population 
asymmetry 

A+- TC B s ,B s -K*-K+i S +0 >S ~o\ ~ TC B s ,B s ^K*+K-l S +0 <S -o\ 0/ , , 

KKtt ~ f s >s ]+T eS l s <s 1 - + 15/ ° • ( 98 ) 

1 B 3 ,B 3 -^K*-K+ L*+0^*-oJ ^ 1 B 3 ,B 3 ~*K*+K- L*+0^*-oJ 

This shows, once again, that focusing on the population asymmetry in the resonance regions on 
the Dalitz plot is sufficient to capture the CP-violating effect of interest. We note that our 
result is close to the value of A CP for B° s — > K* + K~ decay in Table II, and this is caused by the 
dominance of the B° s — > K* + K~ and B® — > K*~ K + contributions to the K*^ bands in the Dalitz 
plot. 



V. CONCLUSIONS 



The study of direct CP violation in the I?-meson system is needed to clarify the nature and origin 
of CP violation in nature. In this paper, we have studied a new tool: the population asymmetry 
in the untagged B d , B s decay rate to a self-conjugate final state of three pseudoscalar mesons. The 
population asymmetry emerges from the failure of mirror symmetry across the Dalitz plot; it signals 
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the presence of direct CP violation. The amplitudes which give rise to such an asymmetry are of 
two, distinct classes. The first, proposed in Ref. [23], follows if we can separate the self-conjugate 
final state, via the resonances which appear, into distinct, CP-conjugate states; we term the latter 
CP-enantiomers. We note that such states are distinguishable by their location with respect to 
the mirror line of the Dalitz plot; in this sense they are distinguishable, mirror images. In this 
construction the angular momentum of the resonance with respect to the bachelor meson is fixed; 
such amplitudes form part of a more general class of amplitudes, for which the angular momentum of 
the two interfering amplitudes are of the same parity, as delineated in Eq. (42). The second class of 
amplitudes emerges from interfering amplitudes whose angular momentum is of differing parity, as 
delineated in Eq. (43). As an explicit example of the latter, we studied the CP- violation consequent 
to s-wave and p-wave interference, mediated by the p and a resonances, in B — > 7r + 7r~7r° decay. 
In our numerical examples, we find that the population asymmetry in the untagged decay rate is 
typically comparable in size to the partial rate asymmetries, making its measurement statistically 
advantageous. 

The population asymmetry serves as a complementary assay to the study of the time-dependent 
asymmetry in neutral 5-meson decays. In this regard, we have emphasized the SM tests discussed 
in Ref. [53], specifically the equality of the quantity Sj, induced by the interference of B-B mixing 
and direct decay, for different modes mediated by b — > sec, b — > sss, and b — > dec transitions. It 
is complementary in that it can simultaneously be regarded as a test that "polluting" SM ampli- 
tudes are indeed as small as expected, as well as a discriminant between different models of new 
physics, should the deviations from the SM predictions turn out to be large. These ideas find useful 
application in our discussion of B — > K S K + K~ decay. One generally expects that physics beyond 
the standard model is most likely manifested through loop-level effects in B decays [53, 64]; new 
contributions may well enter neutral- B mixing or the penguin part of the decay amplitudes, or both. 
Consequently, in the CP-violating observables we consider, new physics is most likely relegated to 
the penguin terms. In this regard, the determination of Ap~ 5 in B s decay may be the most 
sensitive discriminant of emergent phenomena, as its asymmetry is expected, on general grounds, 
to be very small in the SM. 
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APPENDIX 

For the B — > pit amplitudes in Eq. (1), the relevant matrix elements in the framework of Ref. 
are 

(n + (p)\u % Ld\0) = (n-(p)\d % Lu\0) = if nPfl , 
(p + (p,e)\u % d\0) = (p-(p,e)\d % u\0) = f p M p£ ; , 
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q»(p + (p,e)\u % Lb\B°(k)) = -2iA^(q 2 )M p e*-q } 

M 2 - M 2 (100) 
(n+(p)\uYLb\B°(k)} = (k + prFf^tf) + B ^ * q» (f^ \q 2 ) - F*~** \q 2 )) , 

where L = 1 — 7 5 , q = k — p, f n and / are the usual decay constants, and A (g 2 ) and i*o i(<Z 2 ) 
are form factors. In our phase convention for the meson flavor wave-functions, B° = bd, n + = ud, 
\^2tt° = uu — dd, Ti' = du, , with the same convention for the p. Then, using isospin symmetry, we 
also have (n + \u^Ld\0) = V2 (ir°\u^Lu\0) and (<K + \u^b\B°) = -y/2 (n Q \d^b\B°) . The resulting 
amplitudes for B° — > pn are 

a + _ = V2G F {\W ai -\ 



(d) 



''6 



}LM p A^(M 2 ), 



a_ + = ^2G F [A? a, - A?> (a 4 + a 10 )] f p M p F^{M 2 p ) , 



v/2 
V2 



+ A 



(d) 



a 4 + |a 7 - |a 9 - ±a 10 - (a 6 - ±a 8 ) i^o] } f w M p <^ P (M 2 ) 
|a 10 )]/ p M p Ff^(M p 2 ) , 



(101) 



2 a 7 2 °9 



where the numerical values of a i for 6 — > d decay are obtained from Table III, i?^- = 2M 2 _/[(m M + 
m b )(m u + ra d )], and -R^o = M^o/[(w d + m 6 )m d ]. The amplitudes a fei for £>° — > p7r can be found 



from a ife by replacing A^ with A^ J * and using the b — > <i entries in Table III. 

Turning to the £> — > p7r — > 37r amplitudes in Eqs. (2, 4), 9 we note that they contain the p — > 7r7r 
vertex function r p7r7r (s), which is given by 



,(<?)* 



pirn 



/ (TV 



(102) 



where F (s) is the vector form-factor of the pion and / is the p-j coupling constant. The form 
factor F (s) is determined by fitting to e + e~ — > 7r + 7r~ data with a parametrization consistent with 
theoretical constraints [65]. The value of / is determined from the p — > e + e~ width, which, in 
turn, is extracted from the e + e~ — > 7r + 7r~ cross section at s = M 2 [66]. The overall sign is chosen 
so that Eq. (102) is consistent with the Breit-Wigner (BW) form, 



■pBW 

pTTTV 



1 p7T7V 



s - M 2 + iT p M p 



(103) 



as s -> M 2 , where g p7Tn is the p — > 7r7r coupling constant. In our numerical analysis, we employ the 



"solution B" fit of Ref. [65] for F p (s), for which f py = 0.122±0.001 GeV 2 [66]. For s > (M n + M UJ ) 



The signs of the different terms in Eq. (4) follow from the p — > -kit couplings (71-° (p ) 7r ± (p ± )|p ± ) = ±g p7T7r e p ■ 
(Pq ~ P±) an( i ( n+ {p+) n ~{P-) \p°) = 9pn-K e P ' (p+ ~ where we adopt the notation (Jl^il^lMj) = 

^M 2 M3|7i strong |M 1 ^. The signs of these couplings follow, in turn, from the phase convention we have chosen 
for the flavor wave- functions: |7r ± ) = =p| J = 1, 1 3 = ±l) and \ir°) = [/ = 1, J 3 = 0), etc. 
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TABLE III: Numerical values of the coefficients a i with iV c = 3, reproduced from Tables VI and VII in 
Ref. [38], for b — > q and b — > q transitions, with q = d,s. Each of the entries for a 3 , • • • , a 10 ought to be 
multiplied by a factor of 10~ 4 . 





b^d 


b^d 


b^s 


b^s 




1.05 


1.05 


1.05 


1.05 


a 2 


0.053 


0.053 


0.053 


0.053 


a 3 


48 


48 


48 


48 


a 4 


-412- 36 i 


-461 - 124 i 


-439 - 77 i 


-431 - 77 i 


a 5 


-45 


-45 


-45 


-45 


a 6 


-548 - 36 i 


-597- 124 i 


-575- 77 i 


-568 - 77 i 


a 7 


0.7- l.Oi 


0.3- 1.8 i 


0.5- 1.3 i 


0.5 1.3 i 


a 8 


4.7- 0.3 i 


4.5- 0.6 i 


4.6- 0.4 i 


4.6 - 0.4 i 


o 9 


-94- l.Oi 


-95- 1.8 i 


-94- 1.3 i 


-94- 1.3 i 


a io 


-14- 0.3 i 


-14- 0.6 i 


-14-0.4i 


-14- 0.4i 



we match T pnjT (s) to the BW expression, as the fit for F p (s) exists exclusively in the region where 
7T7T scattering is elastic, see Refs. [51, 52] for further details. 

To incorporate the B — > an contributions, we need the additional matrix element 

q»(a(p)\d % Lb\B°(k)) = -i [M\ - M 2 ) f B ^(q 2 ) . (104) 
It follows that the B° — > an amplitude in Eq. (67) is 

Gr, 



{A? a 2 + Af° [a A + |a 7 - |a 9 - \a w - (a 6 - |a 8 ) it^o] } (M^ - M 2 ) /^'(M 2 ) 
- G F (a 6 - |a 8 ) (Ml„ - M 2 ) F^M*) , (105) 



where we use M a = 478 MeV determined in Ref. [67], F B "^(M 2 ) = 0.46 as per Refs. [50], 
and (a\dd\0) = 2M 2 + /[ v / 6 x (m u + m d )}, with X = 20.0 GeV" 1 [52]. To obtain the B° -> (jtt 
amplitude, we note that CP\p°TT°) = + |p°7r°) and CP|cr7r°) = — |o"7r°) for B — > p7T,o-K, as the 
relative angular- momenta in the final states are Z = 1 and 0, respectively. Consequently, a CT can 
be found from a CT by replacing A^ with A^*, using the b — > d entries in Table III, and adding an 
overall minus sign, which is consistent with Eq. (71). 

The a contribution to the B — > 3n amplitudes in Eq. (68) involve the function f a (s) = T an7T (s), 
which describes the a — > nn vertex function in the isospin 1 = and angular-momentum J = 
channel. This function must satisfy known theoretical constraints [52], and we adopt 

r_(*) = xTT(s) , (106) 

for our numerical work, where the function T"(s) is computed in Ref. [68]. We remark that the 
theoretically consistent T UJTn (s) is very different from the Breit-Wigner form used in Ref. [67], as 
detailed in Ref. [52]. 
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For the B — > D* ± D =F amplitudes in Eq. (80), the needed matrix elements are analogous to 
those of the /wr case, where the replacement of the u with c quarks is realized throughout. Thus for 
B° — > D*D decay, we have 



a + _ = V2G F {A«> a, - \f [a 4 + a 10 - (a 6 + a 8 ) i? D ] } f D M D , A^ D * (M 2 D ) 
a_ + = V2G F [A^ fll - Aj d) (a 4 + a 10 )l / D ,M D , F^ D {M 2 D A , 



(107) 



where i? D = 2M|,/[(m c + m b )(m c + mj. For the form factors, we use [69] 



2F^ D {q 2 ) = " >n ' '"'> fjw) , 2A^ D ' 1 ^ - m B + m D* 



{<?) 



B D 



m n m 



(108) 



B D* 



where 



1 + w 



w 



m\ + m 2 DM - q 2 
2m B m DW 



(109) 



The amplitudes a kl for B° — > decay are obtained from a ;fc by replacing A^ with A^* and 

using the b — >■ d entries in Table III. 

The B — > D ± D Tr n° amplitudes in Eq. (81) contain the functions 



/- 



^2 9d*Dtt 



± 



s + _ - s T0 -m 2 D + ml- (m| - s ±0 - m 2 D ) 



m- D -m; 



m 



(110) 



where g D * Dw is the coupling constant associated with the strong decays D*^ — > D^n , 



(111) 



From the empirical values £(£>*+ -> D+tt ) = (30.7±0.5)% [70] and = (96±4±22)keV, [71] 
we extract g D * Dw = 6.3. The use of the Breit-Wigner form for the D* resonance in Eq. (110) is 
appropriate as the D* is narrow. 

The B s — > K^K^ amplitudes also follow from the pir case, and thus we have for B° s — * K^K^ 



V2G F {a« fll - A« [a 4 + a 10 - (a 6 + a 8 ) i^] } f K M K , A^ K \M 2 K ) , 



V^G F Ai s) a x - \\ s) (a 



l 10y 



(112) 



f K ,M K ,F^ K (M 2 K ,) , 



where the values for b — > s decay are given in Table III and _R X = 2M| r /[(m u + m fe ) (m M + m s )]. 



The amplitudes a kl for B° — > K* ± K T follow from a lk by replacing A^r with A^, and using the 
b — > s entries in Table III. The K* ± — > fT 1 * 1 ?]- coupling constant occurring in the analogue of f ± in 
Eq. (110) is g K * KlT = 3.2, which is extracted using the measured K* + width [70] and the isospin 
relation T Rt+ ^ Kn = 3T Kt+ ^ K+7r0 . 
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Finally, we list all the input parameters which have not been previously mentioned. For the 
heavy mesons, we have 

M B o = 5.2794 GeV , r B „ = 1.548 x 10~ 12 s , 
Moo = 5.3696 GeV , r B „ = 1.493 x l(T 12 s , 

(113) 

f D = 200 MeV , M D± = 1869.3 MeV , f D « = 230 MeV , M D *± = 2010.0 MeV , 
f Ds = 240 MeV , M D ± = 1968.6 MeV , f D , = 275 MeV , M D » ± = 2112.4 MeV , 

and for the light mesons 



f w = 134 MeV , M n± = 139.57 MeV , M w0 = 134.98 MeV , 
f p = 210 MeV , M p = 769.3 MeV , T p = 150 MeV , g pn7T = 5.8 
f K = 158 MeV , M K± = 493.68 MeV , 
f K . = 214 MeV , M K *± = 891.66 MeV , T K , ± = 50.8 MeV , 



(114) 



where the decay constants have been obtained from Ref. [72]. For current quark masses, we use the 
running values at /i — 2.5 GeV found in Ref. 



m u = 4.2 MeV , m d = 7.6 MeV , m s = 122 MeV , 

(115) 

m c = 1.5 GeV , m b = 4.88 GeV , 

For the heavy-to-light form factors, we employ 

A^ p {Ml) ~ A^ p (0) = 0.28 , F^{M 2 p ) ~ F^(0) = 0.33 , 

F^{M 2 a ) ~ F B ^(0) = 0.33 , (116) 
A*^ K *(M 2 K ) ~ A*^ K \0) = 0.24 , F*^ K (M 2 K *) ~ F?'^ K (0) = 0.27 , 

estimated in the BSW model [73]. 
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